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Abstract 

Various integrable geodesic flows on Lie groups are shown to arise by 
taking moments of a geodesic Vlasov equation on tlie group of canon- 
ical transformations. This was already known for both the one- and 
two-component Camassa-Holm systems |GiHoTr200"5] IGiHoTr2007l . The 
present paper extends our earlier work to recover another integrable sys- 
tem of ODE's that was recently introduced by Bloch and Iserles [Blls2006] . 
Solutions of the Bloch-Iserles system are found to arise from the Klimon- 
tovich solution of the geodesic Vlasov equation. These solutions are shown 
to form one of the legs of a dual pair of momentum maps. The Lie-Poisson 
structures for the dynamics of truncated moment hierarchies are also pre- 
sented in this context. 



1 Introduction 

Kinetic equations govern the evolution of probability distributions in the phase 
space of many-particle systems in non-equilibrium statistical mechanics. For ex- 
ample, the phase-space probability distribution of a many-particle system whose 
correlations are negligible is governed by the collisionless Boltzmann equation, 
also known as Vlasov equation |V11961| . This equation encodes evolution of the 
Vlasov single-particle probability distribution /(q, p, t) as conservation along 
phase space trajectories, written as 

df ^ , . dh . dh 

where h is the single-particle Hamiltonian, often expressed as the sum of kinetic 
and potential energies 

/,^M! + T/(q) 

in physical applications. Applying the canonical Poisson bracket {•, •} in the 
phase space variables (q, p) expresses the Vlasov equation in its familiar form 

f + {/-/^}=0. (2) 

The Hamiltonian structure of this system is well known |MaWel981| . Namely, 
the Vlasov equation possesses a Lie-Poisson bracket defined on the Lie algebra 
of canonical transformations, such that the Liouville theorem for preservation 
of the volume on phase space entirely characterizes Vlasov dynamics. 

The moment method, widely used in probability theory, provides approxi- 
mate descriptions of the Vlasov solutions. Moments are functionals of the distri- 
bution function / obtained by projections onto the space of phase space polyno- 
mials (symmetric tensors) . Since the Vlasov distribution depends on both posi- 
tion q and momentum p, one may define two different types of moments. These 
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are the kinetic moments and the statistical moments. Kinetic moments are given 
by projection of the Vlasov distribution /(q, p) onto the n-th power of the single- 
particle momentum, p. In contrast, statistical moments are integrals of /(q, p) 
against the n-th power of the full phase space vector, z = (q, p). The remarkable 
property of these two hierarchies of moment projections is that they each define 
equivariant momentum maps |HoLyScl9"90l IScWel994"l IGiHoTr20'08] . Conse- 
quently, the resulting moment dynamics is again Lie-Poisson. Moment equa- 
tions possess interesting closures, which are given by the particular Lie algebra 
structure determining their Lie-Poisson bracket. For example, kinetic moments 
of the Vlasov equation at zero-th and first order yield the familiar closure known 
as ideal fluid dynamics. 

Remarkably, these kinetic moment equations are associated to a family of 
integrable dynamical systems, whose most famous example is probably the Ben- 
ney system for shallow water dynamics jBel973| IGil981j . This convergence of 
different areas of mathematical physics also occurs for several other integrable 
equations. For example, as shown in [GiHoTr2005| . a specific form of the first- 
order kinetic moment equation yields the Camassa-Holm equation |CaHol993] . 
Extending the system to include the zero-th order moment yields another inte- 
grable system; the two-component Camassa-Holm system |CiHoTr2007| . Inter- 
estingly enough, these Camassa-Holm systems are geodesic flows on different Lie 
groups, arising as moment closures of the same kinetic equation, called here the 
geodesic Vlasov equation or EPCan. The latter acronym refers to the Euler- 
Poincare (EP) equation on the group of canonical transformations Can acting 
on phase space r*(50 A special case of EPCan for canonical transformations 
whose generating functions are linear in the canonical momentum has recently 
appeared in the theory of metamorphoses in imaging science |HoTrYo2007| . 

This paper reviews the theory of geodesic equations on the statistical mo- 
ments and shows how such equations possess an additional interesting closure, 
which is related to the space of purely quadratic Hamiltonian functions. We find 
that such a closure yields a particular case of yet another integrable system, re- 
cently discovered by Bloch and Iserles |BlIsMaRa2005[ IBlIs2006] . Moreover, ex- 
tending to inhomogeneous quadratic Hamiltonians yields complete equivalence 
between moment equations and the Bloch-Iserles (BI) system. 

Plan The rest of this section adds a few more remarks about the Lie-Poisson 
bracket for the Vlasov equation. Section 2 is devoted to the Hamiltonian struc- 
ture of the Vlasov moments and their truncations. Section 3 formulates the 
geodesic Vlasov equation, presents its dual pair and illustrates the geometric 
footing of kinetic theory. Section 4 shows how both of the Camassa-Holm sys- 
tems are obtained as geodesic equations on kinetic moments. The last section 
derives the BI system from the statistical moment equations and presents the 
corresponding Klimontovich solutions. 

^The group of canonical transformations Can is also known as the Hamiltonian diffcomor- 
phisms Diffn^^. 
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1.1 The Vlasov kinetic equation 

The Vlasov equation is a Lie-Poisson Hamiltonian system on the group of canon- 
ical transformations of the phase space T*Q for a configuration manifold Q 
|MaWel981| . The dynamics of Lie-Poisson systems takes place on the dual g* 
of the Lie algebra q of the symmetry group G. In this case G = Can{T*Q) and 
= ^can{T*Q). That is, the Lie algebra is the infinite-dimensional space of 
Hamiltonian vector fields. Given the Lie algebra isomorphism between Hamil- 
tonian vector fields and phase-space functions (Xcan — J'), the dynamical 
variable is a phase-space distribution /(q, p), i.e., a density on phase space 
(/ € ~ Den). Upon using the definition of canonical Poisson bracket {•, •}, 
the Vlasov Lie-Poisson structure is found to be 

{F^H}[f]^ ||/(q,p)|^,^|d^'qd^p (3) 

where K = dim(Q). The Vlasov equation ^ is recovered upon choosing F = f 
and h = SH/6f. 

In many physical applications, the Vlasov Hamiltonian is the sum of kinetic 
and potential energy. For example, electrostatic or gravitational interactions 
are governed in the absence of collisions by the Poisson- Vlasov system whose 
Hamiltonian is given by 




where A ^ denotes convolution with the Green's function of the Laplace oper- 
ator. 

2 Hamiltonian structure of Vlasov moments 

The moment method is a popular approach in kinetic systems theory. This 
approach is justified geometrically because taking moments of the Vlasov dis- 
tribution is a momentum map |HoLyScl990l [ScWcl994l IGiHoTr2008) . This 
momentum map arises via the dual of a Lie algebra homomorphism arising 
from the well-known isomorphism between symmetric tensors and polynomials. 
The main point is that this momentum map endows the space of symmetric ten- 
sors with a Lie bracket, thereby generating a well defined Lie algebra. In what 
follows, we shall analyze the cases of kinetic and statistical moments separately 
and then discuss their similarities. 
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2.1 Kinetic moments and the Schouten concomitant 

Kinetic moments are constructed from the following fiber integral |QiTa2004] 

v4„(q, t) / (p • dq)" /(q, p, t) d^'q A d^p 

p,,...p,^ dg*i ®---®dg'"/(q,p,t) d^qAd^p 



ii...i„ = l "'TqQ 



K 



J2 {M<l:t)).,„„,^ dg'^®. 

2l . . An—l 



idg*"Od^q, (5) 



where p • dq denotes the canonical one form (canonical momentum) and d^'q is 
the volume element on the configuration space Q. This construction projects the 
Vlasov distribution onto the space of symmetric tensors. In particular, kinetic 
moments are defined as symmetric covariant tensor fields carrying the volume 
element. That is, they arc symmetric covariant tensor densities. 

The moments arc functionals of the Vlasov density /. Hence, their varia- 
tional derivative may be computed by applying the chain rule as 



5/ ^ ^ 6An 5f ^ ^ 



5{,An)ii...in 



5F 



K 



71 — ii...in- 



SF 

S{An)ii...i„ 



^ SF 



n=0 

which explicitly defines the contraction operation _l . This chain rule formula 
expresses the Lie algebra homomorphism (isomorphism) from symmetric tensors 
to polynomials, whose dual is the momentum map associated to the moments 
|GiHoTr20 08. . Inserting the chain rule formula into the Vlasov bracket yields 
the Lie-Poisson bracket for moments. 



{F,G} [A] 



in which the bracke10 
SF SG 



5An ' 6A„ 



'E 

n.m—O 



5F 



A 



(q) -I 



5F SG 



SA„ ' SAr, 



d^q 



5A„ 



• V ® 



(5G 
SA„, 



5G 

SA^r 



■ V 



5F 

Ja^ 



(6) 



(7) 



is inherited from the canonical Poisson bracket. Here the notation A ■ B for 
one-index contraction between covariant and contravariant tensors is written as 



The operator S here takes the symmetric part of its argument. 
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{A-B)'^j- = Aij,„kB'''^^- and analogously for A = {B-A)'^]^ - = B''"'-^Aiji,„. 
This bracket is well known in differential geometry as an invariant differential 
operator of first order |Nil955j . In fact, this operation is a Lie bracket, which is 
known as the Schouten concomitant or symmetric Schouten bracket. See, e.g., 
|GiHoTr2008| for more discussions and references. 

Remark 2.1 (History of Lie-Poisson structure for kinetic moments) 

In one dimension, the moment Lie-Poisson structure 0j is the Kupershmidt- 
Manin bracket \KuMal978j^ which was found in the context of the integrable 
Benney system in shallow water theory. Lebedev was the first to establish its re- 
lation with the Lie algebra of Hamiltonian vector fields in \Lel979^ and Gibbons 
recognized later \Gil981}j its direct relation to the Vlasov flow. In higher dimen- 
sions, Kupershmidt introduced a multi-index notation \Kul98'fl , corresponding 
to the occupation number representation of the symmetric Schouten bracket. 
This observation suggested the quantum-like framework for kinetic moments in 
]GiIIo Tr2008f . where the moment space is described in terms of a bosonic Fock 
space. 

The moment algebra comprises symmetric contravariant tensor fields and 
these may be characterized as the Fock space represented by a direct sum of 
symmetric powers of vector fields given by 

0-0 V^(Q) ^ith y X:= Si0x\^:3,,. (8) 

71=0 \i=0 I 1=0 \j=0 / 

This is reminiscent of the universal enveloping algebra of the diffeomorphism 
group Diff((5), which is the enveloping algebra of vector fields X(Q) on the 
configuration space Q. It is a standard result that the graded structure of an 
enveloping algebra possesses a Poisson bracket structure [DaSWel999| . 

Remark 2.2 (Kinetic moments and Poisson-Lie groups) 

Interestingly enough, the Schouten concomitant identifies the kinetic moment 
algebra with the Lie algebra of symbols of differential operators. This identi- 
fication is quite suggestive, since differential symbols are known to be a subal- 
gebra of the Lie algebra of pseudo- differential symbols \KhZal 995f . While the 
differential symbols (the moment algebra) are supposed to have no underlying 
Lie group structure, the group ^'D of pseudo-differential operators is a well de- 
fined Poisson-Lie group \KhZal995f . This suggests that the characterization of 
coadjoint orbits for moment dynamics requires the complete Poisson-Lie group 
structure of pseudo-differential operators. A similar direction involving vector 
fields was followed by Ovsienko and Roger \OvRol999f . Also, the appearance 
of the Wick-ordered product from quantum theory in this Poisson-Lie group 
context \KhZal995[ \OvRol999f impl ies a further relation to moment dynam- 
ics, whose quantum-like creation and annihilation operators were presented in 
iGiHoTr2008l 



6 



The moment algebra carries a graded structure (g = (Bi Qi) with filtration 



n+m— 1 (9) 

thereby recovering the space of vector fields gi = X as a particular subalgebra. 
The largest subalgebra is however go® 0i ^X(§)J-, i.e. the semidirect product 
of vector fields with scalar functions. This space appears in the description 
of ideal compressible fluids, where the X- variable is the fluid velocity and the 
JT-variable is associated to the fluid density. Geodesic motion on this space 
has also recently appeared in the metamorphosis process in imaging science 
|HoTrYo2007| . 

The moment equations may be written in Lic-Poisson form as 
dA. 

^TT" = - a,d*_siL An+m.-l (10) 

Ot <5A„ 
??i,n— 

where ad* is the Lie algebraic coadjoint operation defined using the pairing 

oo oo 
/c,n— A:,n— 

The explicit expression for ad* is given in [GiHoTr2008| in any number of di- 
mensions. Here we present the one-dimensional case which will be needed in 
the following sections. The Schouten concomitant in equation ([7]) assumes a 
particularly simple form in ID 

[am,/3n] = ma^dgPn - n(3ndqam (11) 

where q is the spatial variable. Simple use of integration by parts yields the 
following covariant tensor density of rank k — n + 1: 

adl^Ak = (fc + 1) Ak dqPn + n /?„ dqAu . (12) 

This operation was introduced by Kirillov |Kil982| . who first envisioned the 
possibility of a Lie-Poisson bracket on the symmetric Schouten algebra. Familiar 
versions of this operator with n — 1 or fc = n arise in the theory of ideal fluid 
dynamics, soliton dynamics and image matching, while very little is known 
for other values on n, k. Some features of this intriguing open question are 
investigated further below, in dealing with truncations of moment hierarchies. 



2.2 Statistical moments and their Lie-Poisson structure 

As we have seen, the fiber integral defining the kinetic moment hierarchy in 
([5]) requires a kinetic equation on a cotangent bundle. In contrast, the notion 
of statistical moments is given on a symplectic vector space. Upon denoting 
z = (q, p) , the definition of the n-th statistical moment is given by 

X"(t) := / z"/(z,i) d^z (13) 
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where the upper index n in the integrand denotes tensor power (z" = d2)"z), 
while for the time-dependent tensors it denotes the tensor rank. This 

definition places the statistical moments and kinetic moments into the same 
mathematical framework. The first observation is that statistical moments are 
symmetric contravariant tensors on phase space, which is now a symplectic 
vector space V of even dimension N ~ 2K (eventually V ~ M.^) with elements 
z = z'e, e V. 

The moment Poisson bracket for statistical moments may be obtained by fol- 
lowing exactly the same steps as in the previous discussion for kinetic moments. 
That is, one inserts the chain rule formula 



SF 



E 

oo 

E 



(dF\ 








J 


ii . ..in 


V Sf 




dF 


SX" 


OO 

-E 

n=0 


OF 




Sf 





_l z" 



(14) 



(and definition of _l ) into the Vlasov structure, which may then be written as 

d SG" 



{F,G}[f] 



d^z 



f_l — — 

^dz Sf ^ dz Sf 



d^z 



where J is a non-degenerate two form. That is, J is a iV x antisymmetric 
matrix of maximal rank. Since equation (|14p implies 



d SF 
diJf 



E 



dF 



it follows that the moment Poisson structure is 



_l z" 



(15) 



{F,G}{X) 



^ ^ j^n-\-m — 2 
n,m— 



OF dG 



^ E 

n,m— 



X 



n+rn — 2 



■ dF dG 



QXn ' QXr 



where 



dF dG 



QXn ' QX-- 



s 



dF 
dX^ 



dG 
dJc^ 



(16) 



(17) 



is the moment Lie bracket, in which again S operates to take the symmetric part 
of its argument. Recall that J is considered as a contravariant antisymmetric 
matrix, i.e. it possesses upper indexes J*-' = —J-'*. 

Thus, again, the isomorphism between symmetric tensors and polynomi- 
als produces the momentum map associated with the moments jScWel994] . In 
turn, this means that the Lie-Poisson bracket for statistical moments is inherited 
from the Vlasov Lie-Poisson structure. In contrast to the Schouten concomi- 
tant ([7]) for kinetic moments, the Lie bracket for statistical moments in (fTT)) still 
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involves the symplectic matrix J (Poisson tensor). Thus, the dynamics of the 
statistical moments depends explicitly on the original symplectic structure. This 
allows, for example, the direct construction of moment invariants (Casimirs) as 
presented in |HoLyScl990| . Also, the moment algebra ([T7]) involves symmet- 
ric tensors that are covariant, rather than contravariant as happens for the 
Schouten concomitant. 

Remark 2.3 (Statistical moments and accelerator beam optics) 

Statistical moments are important, for example, in the study of beam dynamics 
in particle accelerators. In this framework, they are defined as \ Chi 98 3^ 



M^{t) -.= 11 p"q"/(q,p,t) d^q d^q 

= (-M^(i)) 



^1 



where e^ is a basis element of the configuration vector space Q, while e*^ is its 
dual (so that V ~ Q x Q* ). In ID, the beam emittance 

e:= {MlMt-{Mlff\ 

known as the Courant- Snyder invariant \CoSnl958}/ . was recognized as a mo- 
ment invariant ( Casimir). This observation led to the study of more general 
moment invariants J Dr 19901 \IIoLyScl990 A geometric investigation of statis- 



tical moments was carried out in \ScWel99^^ , where the moment algebra was 
related to the Heisenberg algebra on phase space. For particle accelerator design, 
moments are often used in computational efforts to account for space charge ef- 
fects and other beam-related problems. 

As for the kinematic moments, one characterizes the Lie algebra of statistical 
moments by using the grading, 

oo / n \ oo 

0-0 (18) 

n=0 \ j=0 / i=0 

with the filtration 

[fln ,0m] ^ n+m — 2 

(19) 

which shows how symmetric matrices 02 = ^* V V* = Sym*(A^) form a par- 
ticular subalgebra (here we denote by Sym* {N) covariant symmetric matrices) . 
The largest subalgebra is given by go©0i©02 = M © F*© Sym*(iV) and it 
will play a central role in the remainder of this paper. 

Remark 2.4 (Occupation number representation) 

Analogously to the Lie algebra of kinetic moments, the statistical moments 
also carry a bosonic Fock space structure, where appropriate occupation num- 
bers may be defined through the introduction of the multi-index notation z'^ := 
(z^y^ . . .{z^y . Indeed, taking moments by {t) = J z'^ f{z,t) d^z yields 
the occupation number representation for statistical Vlasov moments. 
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The moment equations are written in terms of the Lie algebra coadjoint 
operator ad* as 



dt ^ — ' ax 

n=0 



=— f:«'5((^.j)jx"— ) (20) 



so that 



n=0 \^ /ji...j„-ik ) 

in which repeated tcnsorial indexes are summed. Having established the nota- 
tion for the various operations among moments, we turn next to their application 
in geodesic Vlasov flows. 

2.3 Truncation of moment hierarchies 

The truncation of Vlasov moment hierarchies is a typical problem in kinetic 
theory and, for statistical moments, this question was addressed by Channell 
|Chl995| by using Levi's decomposition theorem. ChanncU's result may be 
summarized by saying that, for a moment Hamiltonian not depending on the 
first-order moment {dH/dX^ = 0), the moment hierarchy can be truncated at 
any order, thereby yielding a truncated Lie-Poisson system. In order to see how 
this works in practice, let us write the truncated equations at the order K for the 
hierarchy of statistical moments X" = ad,*^ ^ where hm = dH/dX"^. 

We have 



X' = ad;^ X' + ad;3 X^+--- + adl^ x""-' 

X^ = adl X^ + adl X^+--- + adl^ X^ 

X' = adl^X^+ ■■■ + adl^_^X^ (21) 



X^-i = ad;,^ X^-^ + ad,,3 X^ 
X'"^ = ad;;^ X^ 

and we recognize that the equation for X^ is decoupled, so we restrict to the 
equations for X^ , . . . ,X^ . At this point, one verifies that the Lie Poisson 
bracket for the truncated moment system is given by 



n=2 m=2 ^ ^ 



(22) 
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with the same notation as in PT|) . We recognize that the truncated structure is 
completely determined by the Lie algebra filtration p9)) and does not depend 
on the particular expression of the Lie bracket, which was not used in deriving 
the truncated system above. 

Following similar arguments, one can write the truncated system for kinetic 
moments as 

dtAo = -ad,t, Ao - ad^^ Ai ad^^, Ak-i 

dtAi = -ad,*, Ai - ad^^ A2- ad^^, Ak 

dtA2 = -ad;, A2 ad;^„, Ak (23) 



dtAK-i = -ad,;, Ak-^i - ad;^ Ak 
OiAk = -ad;, Ak 

where hm = SH/5Am and we have assumed ho = 0. The zero-th moment 
equation decouples and one is left with the remaining equations for Ai, . . . , Ak- 
These equations possess the following bracket structure: 



{F,G}{A) = -E E 

n— 1 m— 1 



5 An ' M„ 



(24) 



with the same notation as in ([7]). One recognizes again that the truncated 
structure is uniquely determined by the Lie algebra filtration ([9]). This fact 
suggests that a similar approach would also apply to the BBGKY moments of 
the Liouvillc equation (reduced probability distributions) , whose corresponding 
Lie algebra is known to possess a similar filtration |MaMoWe84| . However, we 
leave the investigation of the BBGKY moment hierarchy for another time. 



3 The geodesic Vlasov equation 

We have seen that moment hierarchies are equivalent descriptions of the Vlasov 
equation, which allow for geometric closures of the kinetic system (e.g. the ideal 
fluid closure for kinetic moments). The Vlasov equation is a Lie-Poisson equa- 
tion on the Lie algebra of the group Can(T*Q) of canonical transformations, and 
this property is reflected in the Lie-Poisson structure of moment dynamics. It is 
well known that physical systems with interesting geometric behavior are often 
geodesic flows on Lie groups with respect to a metric provided by the system's 
kinetic energy. The most familiar example is probably rigid body motion, which 
is governed by geodesic motion on SO (3). Likewise, Euler's equations for ideal 
fluids may be interpreted as geodesic motion on the volume-preserving diffeo- 
morphisms Diffvoi(IR^) of the 3D flow domain [Arl966j . Another interesting 
example of geodesic motion is provided by the EPDifF equation [HoMaRal998| , 
which governs geodesic motion on the full diffeomorphism group Diff(M"). In 
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many cases, geodesic flows on Lie groups turn out to be completely intcgrable 
Hamiltonian systems. For example, in a one-dimensional flow domain R, EPDiff 
recovers the Camassa-Holm equation for shallow water waves |CaHol993] . 

The problem of geodesic flow on the symplcctic group {symplecto-hydrody- 
namics) was introduced by Arnold and Khcsin in |ArKel998] . The present 
paper pursues this idea by considering geodesic flow on the canonical trans- 
formations within the context of Vlasov dynamics. In particular, we develop 
a geodesic Vlasov equation called EPCan {Euler-Poincare equation on the 
canonical transformations) as an extension of previous work in [GiHoTr2005| . 

Remark 3.1 (Symplectomorphisms vs canonical transformations) 

The name EPCan refers to the Euler-Poincare equation for geodesic motion on 
the subgroup of the symplectic transformations arising from Hamiltonian vec- 
tor fields. The subgroup Can (which could equally well be called Diffcanj may 
be identified as the group of smooth invertible canonical transformations with 
smooth inverses. These transformations coincide with symplectic transforma- 
tions in simple domains such as R^^. In those simpler cases, the geodesic 
Vlasov equation is also known as FiPSymp ]GiHoTr2005[ \GiHoTr2007l 

As mentioned in the introduction, the idea to investigate EPCan was moti- 
vated by the observation that geodesic equations for kinetic moments were found 
to include the Camassa-Holm equation [GiHoTr2005| . In fact, geodesic moment 
equations arise from EPCan whenever the norm may be expanded as a Taylor 
series. Later, it was recognized jGiHoTr2007j that the geodesic moment equa- 
tions also recover the two-component Camassa-Holm equation jChLiZh2005| . 
The latter is a geodesic flow on a semidirect-product Lie group [Kuzmin2007| . 
In order to explain these issues, we shall introduce the EPCan equation and 
show how it specializes to each integrable case. 

Given a symplectic manifold V of even dimension N = 2K, the EPCan 
Vlasov Hamiltonian is defined by 



where z £ V and the kernel Q is chosen so that it defines an appropriate norm 
on Den(7'). When dealing with moments, wc shall restrict to the special cases 
V = T*Q, with Q a general configuration manifold, and V ~ V, with V a 
symplectic vector space. The geodesic Vlasov equation (aka EPCa7i) is written 
simply as 



which coincides with Euler's vorticity equation in 2D when Q = (—A) ^. 

3.1 Euler-Poincare equations on Hamiltonian vector fields 

In order to understand how the geodesic Vlasov equation arises from an Euler- 
Poincare approach, one starts with an invariant Lagrangian defined on the tan- 
gent space of the canonical transformations C : TCan — > R, which is purely 



m = \ j j f{z)g{z,z') f{z') d^z A^z' = \\\f\\l (25) 




(26) 
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quadratic. By the invariance property, we can write the associated variational 
principle on the Lie algebra of Hamiltonian vector fields as follows 

S f ' L[Xh] dt = (27) 

•I to 

where X/j = JJV/i and 

= ^(^dW {]SQ I \/ h), =: L[h], (28) 

in which Q is taken to be a positive-definite symmetric operator so that L[X.h] 
defines a nondegenerate norm. The Legendre transform 

/ = ^ = div(jQJ V/i) ^ h = (^divJQJv) V (29) 

yields the EPCan Hamiltonian in the Vlasov form 

H[f] = ^{.f,d-\f) (30) 

with 

0~divJQJJV. (31) 

This formula specifies the relation between the geodesic Vlasov equation and 
the geodesic motion on the Hamiltonian vector fields. An interesting case occurs 
when Q is the 'flat' operation QX/j — Xj^ which takes contravariant vectors to 
covariant vectors, so that 

divJ(JJV/i)^ = - A/i. (32) 

Then the operator O reduces to minus the Laplacian 

O = - A (33) 

and in two dimensions one obtains the Eulcr Hamiltonian for vorticity dynamics, 

H[u;] ^ 1/2 {lo, {- A)-' u), 

with to = f. In the more general case when Q is a purely differential operator, 
one finds that Q and J commute and thus O ~ — div Q V. Also if Q commutes 
with the divergence, then, one has O = — QA. However in general, Q is a 
matrix differential operator that does not commute with J. 



13 



3.2 The ideal fluid dual pair for the Vlasov equation 



Both the Poisson- Vlasov system (g]) in plasma physics and the geodesic Vlasov 
equation (|25p possess the single-particle solution (Klimontovich solution) ex- 
pressed as a delta function in phase space, 

^^ 

/(z,t)=^u;,<5(z-CW), (34) 

where the index a is summed over a = 1, . . . , TV and the number Wa is a constant 
weight associated to each particle. This apparently trivial solution is of funda- 
mental importance in physics and leads to the Klimontovich theory of kinetic 
equations |K11982| . In the remainder of this paper, we shall see how such a 
solution emerges within the analysis of integrablc moment closures of EPCan. 

In particular, an interesting situation occurs when one allows for more gen- 
eral solutions of the form 

AT 

/(z, t) = J2 ^-(^) - ^ais, t)) d'^S (35) 

a=l •' 

where s is a coordinate on the immersed submanifold Sa C with k = dim(S') 
(e.g. k = 1 for a curve, or fc = 2 for a surface). For simplicity of notation, we 
have already suppressed the implied subscript a in the arclength s for each Wa 
and . For further simplicity and without loss of generality, we take A/" = 1 and 
so suppress the index a in what follows. This is equivalent to treating an isolated 
singular solution. As one might expect, this is only a notational simplification; 
not a real restriction. 

The C's in ^ belong to the space of embeddings Emb(S',R2^). That is, 
: S ^ M."^-^ . Remarkably, these solutions define a momentum majo 

Jsing : Emb(5, R2^) X^.JR^A) (gg) 

where X*an(R^'^) — Den(R^^'). This momentum map is produced by the left 
action of canonical transformations on V by composition of functions; that is, 

?/ • C = ?? o C • 

Importantly, the same kind of momentum map arises in the motion of ideal fluids 
(e.g., for point vortices in 2D). See |MaWel9"83] . where these momentum maps 
are shown to possess a dual pair structure [Wel983| . Namely, if one considers 
5 as a manifold with volume form 

Wvoi = w{s) d'^s , 

•^Here SEcan denotes the Lie algebra of Hamiltonian vector fields, whicfi should not be 
confused with the larger Lie algebra corresponding to the tangent space at the identity of the 
symploctomorphism group Symp. 
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then the right action of Diffvoi on Emb(S', R^-^) 

C-v = C°v, 

yields another momentum map 

Js : Emb(S', M^^) ^ X^^^iS) . (37) 

In more generaUty, if (S^w) is a volmire manifold and (VyUj) is a symplectic 
manifold, then the right action momentum map is (cf. [MaWel983) ) 

J5 : C ^ C*t^ : Emb(5, V) ^ X*US) . (38) 
To summarize, we have the following dual pair structure 



Emh{S,V) 




which is formally equivalent to the dual pair structure for ideal fluids [MaWel983| . 
Moreover, the left leg yields a solution of the Vlasov equation regardless the 
number of dimensions and this makes the above dual pair a natural object in 
kinetic theory. 

In order to write explicit formulas, we specialize to the case 
denoting C,{s) (Q(s),P(s)), one defines the following Poisson structure on 
Emb(5, V) 

r ^ ^1 ^ f I f SF 5G SG 5F\ , 

= g y - - ^ - j d^. (39) 

where we see that the factor l/w(s) is needed for functionals of the form 
G{C) = / oJvoigiC) = / ^{s) g{C{s)) d''s, whose functional derivative SG/6( = 
w(s) dg/dC takes values in Den(S'). 

Remark 3.2 (Comparison with point vortices) 

The factor l/w(s) is reminiscent of the vortex strength factors in the Poisson 
bracket for point vortices \MaWel983f . Indeed, the bracket above appears as the 
higher dimensional version of the vortex bracket, so that the one dimensional 
vortex strengths are replaced by appropriate densities (the weights 'w(s)) on the 
embedded space S . 

Finally one checks that, for any Hamiltonian fimction h E T{T*R^^), 

{F,{Jsin,,h)}^^^^XH[F] (40) 
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where X/j [F] is the infinitesimal generator of the action of canonical transforma- 
tions Can(M^^) on Emb(S', K^^). Thus, Jsing satisfies the classical definition 
of a momentum map. 

The singular solution momentum map Jsing produces the collective Vlasov 
Hamiltonian 7J o Jgi^g. In particular, substituting the singular solution momen- 
tum map ([55)) into the EPCan Hamiltonian ([^5)) yields the collective Hamilto- 
nian 



,,6=1 



1 r r 

oY. Wais)wbis')g{Qa{s),Pais),Qbis'), Phis')) d" S d" s' (41) 



thereby producing the following collective equations of motion 



dt 5Va 

h—1 



= Wais) ^ / Wb{s') — g{Qais),Pa{s), Qb{s'), Phis')) d'^ s' 
b=l •' 

ais)J2j WUS') ^ g{Qa{s),Pais),Qbis'), Pb(.s')) d" s' . 



dPais,t) ^ _5Hj^ ^ 

dt (5Qa 

b—1 

Remark 3.3 (Possible divergent terms in collective motion) 

The collective dynamics of singular solutions deserves some care, depending on 
the form of the Vlasov Hamiltonian. The existence of such solutions does not 
guarantee the existence of a well defined collective Vlasov Hamiltonian -ffojgjng. 
This is because the singular solution momentum map Jsing may produce diver- 
gent terms in the collective Hamiltonian \Wel983f . For example, this is the 
case of the Vlasov Poisson system Q), where the divergence is generated by 
potential terms such as l/2'^WaWh\Qa — Qb\^^, when a = b. The same sit- 
uation occurs for point vortex solutions of the planar Euler's vorticity equa- 
tion; these solutions correspond to the 2D phase space Hamiltonian Hj\f = 
1/2 '^WaWb log KQa — Qh, Pa ~ ^'h)|- On the other hand, these problems are ab- 
sent, for instance, in the Vlasov-Helmholtz system (see \GiHoTr2008f and refer- 
ences therein), since the potential terms there are given 6j/ 1/2 ^ WaWb e'*^""'^''! . 

As for the right-action momentum map, the expression 

JsiC) = = C*(dq A dp) = Cdq A C*dp = d(C*q) A d(Cp) (42) 

yields the following simple expression 

J5(Q,P) = dQ(,s)AdP(.)= ^ |^^d."Ad.- (43) 

n,m—l 

The conservation law dJs/dt = is recovered by Noether's theorem, due to the 
Diff(5')-invariance of the collective Vlasov Hamiltonian Hj^/ = o Jsing in (|4T|) . 

When S* is a Lagrangian submanifold (this requires dim(S') = l/2dim('P)), 
the momentum map J5 restricts to Js(C) = 0. Likewise, the case dim(5) = 
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recovers the usual Klimontovich solution (|34p of particle motion used in kinetic 
theory. This fact, together with the geometric results on moment hierarchies of 
kinetic equations, illustrates the geometric basis of kinetic theory, in analogy to 
Arnold's formulation of the ideal fluid [Arl966| . 

3.3 Klimontovich solution and the Lagrange-to-Euler map 

This section discusses the two limiting cases of the singular solution momen- 
tum map Jsing, that is dim(5') = and S ~ V. As mentioned above, the first 
case yields the Klimontovich solution (|34[) . which is then a momentum map 
Jging : ^^I'P ^can(^)- ^^^^ casc, thc solution identifies the particle tra- 
jectories, subject to initial conditions Ca(0) = z", so that thc particles are trans- 
ported in the phase space V by canonical transformations as {Ca(0} ~ V't°{za}, 
where tpt S Can{Xa'P) is generated by thc collective Hamiltonian iJ o Jgi^g. In 
two phase-space dimensions, the Klimontovich solution is thc usual point vortex 
solution of thc Euler's vorticity equation. 

Another suggestive case of the above treatment is given by S* = P, so we 
may denote s = z*'. Then, one has Ca(': = "Ht"^ ^ Can{'P) and the momentum 
map Jsing : Can{P) ^can(^) written as 

/(z,t)=^ L,(zO)5(z-,yf'.zO) d^^zO. (44) 

This expression coincides with the well known Lagrange-to-Euler map for fluids, 
whose importance is well established in continuum dynamics. The Lagrange-to- 
Euler map is equivalent to the characteristic form of the Vlasov equation P]). 
Notice that both the Klimontovich and the Lagrange-to-Euler maps are pro- 
duced by the same Lie group Can{V) acting on XaV and XaCan{T') respec- 
tively, with the same left action by composition of functions, that is rj ■ {Ca} = 
{v ° Ca} ™ tlie first case and T] ■ {r/t")} = { 77 o } in the second. On the 
other hand, for thc Klimontovich case, thc collective dynamics generated by 
thc Hamiltonian H o Jgi^g : XaV R- produces thc canonical transformations 
•ip e Can{XaV) 7^ Xg^CaniV). This point is of fundamental importance be- 
cause thc Lie group Can{Xg^'P) is the symmetry group of the Liouville equation 
jMaMoWeM] and the difference between Can{Xg^V) and Xg^Can{V) is related 
to the particle correlations, which are neglected in the second situation. (The 
latter is the Vlasov mean field approximation.) 

The fact that these two fundamental maps each arise from the left leg of a 
dual pair of momentum maps again illuminates the geometric footing of kinetic 
theory. The above arguments also provide mathematical support for the wide 
success of Klimontovich method in kinetic equations |K11982j . 

3.4 Geometric kinetic theory 

The presence of dual pairs in kinetic theory illuminates the Liouville and Vlasov 
equations in thc light of their Lie symmetry properties. (Something similar 
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happens for Euler's vorticity equation.) Namely, the presence of momentum 
maps is not accidental in kinetic approaches. Indeed, a reasonable summary of 
the results in |MaMoWe84l |HoLyScl990l IGiHoTr2008] could be made by saying 
that the process 



Liouville equation — > Vlasov equation 



/ ideal fluid 



beam optics 



(45) 



is given by a composition of momentum maps. In other words, taking the 
moments (BBGKY, kinetic or statistical) of a Lie-Poisson kinetic equation is 
always a momentum map |MaMoWe84l |HoLyScl990l IGiIIoTr2008| . Moreover, 
the closures adopted to obtain Vlasov from BBGKY, fluid theory from kinetic 
moments and beam optics from statistical moments are also momentum maps 
arising from particular subgroups of the symmetry group of the starting sys- 
tem. More explicitly, passing from Liouville to Vlasov requires the subgroup 
Can{'P) C Can{Xi'P). Likewise, passing from Vlasov to fluid requires the fiber 
preserving subgroup CanT^{T*Q) C C'an(T*Q). Finally, passing from Vlasov to 
beam optics requires the subgroup Sp(2/-r, R) C Can(R^^). We can summarize 
the situation in the following statement 



All these moment approximations in kinetic theory are momentum maps. 



That the BBGKY distributions arc momentum maps is a remarkable fact. One 
may ask whether the Klimontovich averages in plasma theory also share this 
property. In this case, the autocorrelations considered in the latter approach 
would again be naturally included in the geometry of the theory. We leave this 
promising question open, as a direction for future research. 



4 Moment closures of EPCan: integrable cases 

As explained in [GiHoTr2005[ IGiHoTr2007j , the geodesic Vlasov equation may 
be represented in terms of the moments. Indeed, upon supposing that the metric 
G in ((25)) and (|26p is sufficiently smooth, may can expand 5(z,z') in a Taylor 
series. 



4.1 Integrable closures of kinetic moments 

In this section, we present the kinetic moment hierarchy for EPCan. Upon 
denoting z = (q, p), one may expand Q in a Taylor series, as follows, 

oo 

g{z,z') = P"®P'"-I G„™(q,q') 

Ti,m— 

= E E E (p").„...,„(p""),„...„jG_(q,q'))--^-— 

n.m 
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where G„„i is now a contravariant tensor field of rank n + m. Inserting this 
Taylor expansion in the EPCan Hamiltonian yields (with the notation above) 

,00 

i?=2 ^ / A„(q)(g.^„(q')-lG„™(q,q')d^qdV=2||Mn}||G (46) 



n,m— 

Thus, upon denoting 



y"G„,„(q,q') -I A„(q') dV (47) 



the geodesic moment equations become 
dA °° 

-Qf^-T. adG^,^,A„+™_i (48) 

m,k—0 

where ad* is the coadjoint Lic-Schoutcn operator. The singular solutions of the 
Vlasov moment hierarchy may be expressed in the following form: 

A„(q, t)^ J w{s) P"(.s, t) d{q - Q{s, t)) d'^S (49) 

which is a momentum map J : Emb(S', T*(3) 0*, where g is the kinetic 
moment algebra. 

As we have seen in the preceding discussions, the moment algebra possesses 
the important subalgebra gi — X{Q) of vector fields on the configuration man- 
ifold. In terms of canonical transformations, this corresponds to Hamiltonian 
generating functions that are linear in the momentum coordinate, i.e. point 
transformations. These are cotangent lifts T*DiS{Q) of diffeomorphisms on 
the configuration manifold Q jHoMa2004| IGiIIoTr20"07| . Remarkably, when the 
moment hierarchy of EPCan is closed such that Gn =: Gi is the only non- 
vanishing term of Gnm, we obtain the Hamiltonian on the one-form density 
Ai =: m(q) ■ dq ® d^q G X* 

H=\jj m(q) • Gi (q, q') m(q') d^^q d^^q' (50) 

By using the property of the Schouten bracket [/3i, a„] = £^3^ a„ (where £ 
denotes Lie derivative), one finds the EPDiff equation, 

+ i^Gi*mm = 0. (51) 

EPDiff is the Euler-Poincare equation on the diffeomorphisms |HoMa2004] . This 
equation has the important property of exhibiting emergent singular (5-like solu- 
tions from any confined smooth initial configuration. In ID, the particular case 
Gi = (1 — q;2(?^)~^ yields the integrable Camassa-Holm equation, which is well 
known in the community of integrable systems. 
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Remarkably, if we also allow for Goo =: Go 7^ 0, we obtain the Hamiltonian 

H[^o,^i] = ^ jjA,{c^)-G,M)A,{c{) d^q dV 

+ ^JJAo{ci) Go(q,q') Ao(q') d^q d^q' 

which yields a geodesic flow on the semidirect-product Lie group Diff(S) J-, intro- 
duced in [HoTrYo2007| in the context of image matching, discussed further in 
|HoTr2009] and shown in numerical simulations to exhibit emergent singularities 
in both of its variables |HoOnTr2009] . Interestingly enough, the EP(Diff(S) J^) 
equations (with notation Pq o Aq = — ad^^ ^0) 

Ao,t + £gi*Ai Aq ~ Q , 

Aij + £g^*a^ Ai = AoO (Go * Ao) . 

are a geodesic flow on the extended point transformations, i.e. compositions of 
cotangent lifts and fiber translations |MaWeRaScSpl983] . 

Remark 4.1 (Geodesic motion on fiber-preserving transformations) 

The semidirect-product Lie group Diff (§)J- is identified with the compositions 
of cotangent lifts with fiber translations on the phase space T*Q with coordi- 
nates (q, p). This identification is of fundamental importance in plasma physics 
\Ma WeRaScSpl983^ . It also yields the interpretation of EP{DiS(§)T) as a 
geodesic motion on the Lie group Can.^{T*Q) of fiber-preserving canonical 
transformations on the cotangent bundle T*Q. In fact, any transformation by 
Can.^{T*Q) can be realized as the composition of a fiber translation and a cotan- 
gent lift (or viceversa) \BaWel997l . On the other hand, such a transformation is 
always a canonical transformation characterized by a generating function that 
is linear (and inhomogeneous) in the canonical momentum. Therefore, since 
Difr(Q)(S) J'(g) Can^{T*Q) then EP(JDiS(§)T) ~ EPCan^. 

In ID, the special case Gq = (5-function yields yet another integrable system. 
At = - (uX)^ 

Ut ~ Uqqt = -iuUq + 2UqUqq + UUqqq - XXq 

known as two-component Camassa-Holm equation [ChLiZh2005| (here u = (1 — 
d^)~^Ai and A = Aq). This system first appeared in |01Rol996] . Its singular 
solutions were studied in |CoIv200"8] who pointed out their relation between 
this system and shallow water equations. Upon slightly modifying the Hamil- 
tonian by Go = (1 — 9^)^^, one has the spontaneous emergence of the Klimon- 
tovich solution ([^5]) in {Aq, Ai), as shown in [HoOnTr2009j . The two-component 
Camassa-Holm system and its Klimontovich solutions have also found applica- 
tions in the metamorphosis approach to image matching, e.g., for magnetic 
resonance images jHoTrYo2007) . 
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4.2 Geodesic flow on statistical moments 

Section 5 shows that representing the geodesic Vlasov equation in terms of 
statistical moments in a special case recovers the well-known Bloch-Iserles inte- 
grable system |BlIsMaRa2005l rBlIs2006l . li G in ^ and ^ is assumed to be 
analytic in both position and momentum, its Taylor expansion may be written 
as: 

oo 

g(z,z')= r™„_lz"®z"" (52) 

n.??i— 

where Tnm is a covariant tensor of rank n-\- m. Then, the EPCan Hamihonian 
(PSj) is written in terms of the moments X'^^{t) = /z^/(z,i) d^z as 

oo 

^=2 ^ r„™_lX"®X" (53) 

n,m— 

so that the moment equations for (n, m) G Z+ tmm out to be 

= -m^nSU (r„fe _l X'') ■ j) _l J (54) 

where, as before, the symbol _l denotes contraction between upper and lower 
indices of the various tensors. In explicit index notation this reads as 

n=0 ^ ^ 

An example of this equation is discussed in Section [5] for the case m — n = 2. 

4.3 Klimontovich dynamics for statistical moments 

As in the case of kinetic moments, the single-particle Klimontovich solution 

/(z,t) = ^i«,5(z-C(i)) (55) 

a 

of the geodesic Vlasov equation offers an interesting opportunity for presenting 
solutions of the moment hierarchy. These solutions may be represented as 

X^t)^J2^aCit) (56) 

a 

where ^" := Ca ® ' ' ' ® Ca times) and Ca satisfies Hamilton's canonical 
equations with Hamiltonian 

a,b n,m 
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This means that satisfies 



n Wa 



b 



Wb Tnm -I Co 



(58) 



where the moment order m or n ranges from zero to infinity, and a, h range over 
the number of particles a, 6 = 1, . . . , A/". Each Co undergoes Hamiltonian dy- 
namics because the Vlasov single-particle solution is an equivariant momentum 
map. Equivariance of the corresponding momentum map for dynamics of point 
vortices in an Euler fluid was proven in |MaWel983| . The same argument holds 
for the singular solutions of kinetic moment equations, including the EPDiff 
equation [HoMa2004] . 

Remark 4.2 (Truncation of geodesic moment hierarchies) 

It is important to notice that higher- order moment truncations do not possess 
the Klimontovich solution. The latter exists only for genuine moment closures, 
such as the fluid closure for the kinetic moments (consisting of Aq, Ai), or the 
2^'^ -order closure for statistical moments used in linear beam optics (consisting 
ofX\X^). 

5 Bloch-Iserles system as a moment equation 

The case (n, m) = (2, 2) of the moment equation (j54p yields an important mo- 
ment subalgebra, given by homogeneous quadratic polynomials on phase 
space, i.e. quadratic forms on V . In this case, formula p7|l implies the following 
Lic-Poisson moment bracket 



{F,G}{X') = A(X\S 



dF , dG 



QX2 - QX2 



(59) 



where the moment X^ is now a N x N symmetric matrix. Because of the 
antisymmetry of J, the bracket above may be rewritten as 



{F,G}{X^) = (X' 



dF dG 



dX^ ' dX^ 



2JI/ 



(60) 



where 



OF dG 



dX^ ' dX^ 



dF _ dG 



dG 



dF 
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QX2 QX2 dX^ dX^ 



(61) 



which is the Lie bracket for the integrable Bloch-Iserles (BI) system of equations 
introduced in |BlIsMaRa2005l IBlIs2006] . In that case, an antisymmetric matrix 
N of any dimension defines the following Lie bracket on the space of symmetric 
matrices of the same dimension. In formulas one has 



[X, Y]n-~ XWf - YNX 



(62) 
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where X and Y are symmetric matrices. The Bloch-Iserles system 

X=[(X)2,N] (63) 

is Lie-Poisson on this Lie algebra, with Hamiltonian Hbi = ^Tr{*XX). Here 
the notation (X)"^ denotes standard matrix muhiphcation of X by itself, in 
order to distinguish from second-order moments X'^. In addition, *X denotes 
the transpose of the matrix X, so that *X ~ X when X is symmetric. One 
concludes the following. 

Theorem 5.1 In the even-dimensional case, the integrahle Bloch-Iserles system 
is the Vlasov moment equation I120\) obtained from the quadratic Hamiltonian 

(64) 



associated with the antisymmetric matrix N := 2J. 

Let us now look at the moment bracket for functions of 

{X°, X\ ^2) e 05 © 0* © 0^ ~ M © V® {V V V) (65) 

where V is again the symmetric tensor product defined in ([5]) and X° = const 
is naturally taken as the probability normalization. The moment bracket (|17p 
becomes 



{F,G}{X)^X'- 



dF 






dX^ 




1 dF 






\dX^ 
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dG 



dG 



dF 



dX^ dX^ dX' 



X' 



dF dG 



dX^ ' dX^ 



(66) 



21/ 



which is given by the direct sum of the canonical Poisson bracket on V in first 
term, plus the semidircct-product Lic-Poisson bracket on 02® 0* ~ Sym(S) V 
in the second and third terms. Thus, the specialization of the moment bracket 
(fT7|) to (|66p in this case defines a Lie-Poisson bracket on (Sym(s)y) ffi R. 

We now turn our attention to the odd-dimensional BI system. In this system, 
one has a degenerate antisymmetric matrix N of odd dimension n and rank 2K. 
Upon defining N as the non-degenerate minor of maximal dimension (2A'), the 
degenerate matrix N 

[N 



N = 







produces the Lie bracket ((62|) associated to the BI equation ([63]). 

The odd-dimensional BI system is known to be a geodesic flow on the space 
{Sym{2K)®M2Kxd) ®Sym{d) endowed with the Lie bracket (cf. equation 
(2.14) in |BlIsMaRa20"05] ) 



[(S,A,B),{S',A',B')] := 
(SNS' 



S'NS, S NA' -S'NA, *A NA' - *A' NA) (67) 
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where one denotes d = n — 2K, for any S, S' G Sym(2if ), A, A' E M-zKxd, and 
B,B' e Sym(d). 

Wc will show that the bracket (|17|) for the Vlasov moment system is Lie- 
Poisson on the dual to the Lie algebra (Sym(2X) (S) y) © R. That is, we choose 
d = 1, which is the case of interest here. Proposition 2.5 in [BlIsMaRa2005) 
shows that the geodesic equations on Sym(2/C + 1) are equivalent to the geodesic 
equations on (Sym(2isr) (S)^) © R, where V is a 2X-dimensional symplectic 
space carrying a non-degenerate symplectic structure N/2. This leads us to 
the identification of the Bloch-Iserles system with geodesic moment dynamics. 

Theorem 5.2 When d ^ I, the odd- dimensional BI system is a Vlasov moment 
equation of the form h2U\) on the statistical moments. This equation is generated 
by the quadratic Hamiltonian 



(68) 



yielding a geodesic moment flow. The corresponding BI Hamiltonian Hbi{X) = 
l/2Tr(*XX) is written in terms of the symmetric {2K +l)-dimensional matrices 
of the form 



X 



X' 



X' 



These matrices are endowed with the Bloch-Iserles Lie bracket \62j). where the 
(degenerate) antisymmetric matrix N takes the form 



N 
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and J is the canonical symplectic matrix. 
Proof. Upon denoting 

dF dG 



N = 2J, 



5 = 



9X2 



S' 



dX^ ' 



A 



dF 



A' 



dG 
dX^ 



(69) 



one sees that the only difference between the Lie bracket in (p7|) and the bracket 
in (l66l) resides in a constant factor in the first term of (l66l) 



dF 



• J- 



dG 



dX^ / dX^ 



dF 
'dX^ 



•2J 



dG 



dX^ \dX^ 



dG 



dF 



dX^ 



(70) 



where the square bracket in the right hand side is identical to the last component 
of (j67p . This difference however can be easily overcome. Indeed, one can always 
re-define the Lie bracket on (Sym(2ii') (S) F) © M as 



[{S,A,B),{S\A',B')] 



sns' -s'ns, snA' -s'nA, - (Una' - *a' na) 
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and verify that the map 



* : (Sym(2is:)(S)l/) ®R 
iS,A,B) 



Sym(2ii: + 1) 

S A' 

*A AB 



is a Lie algebra isomorphism for Sym(2if + 1), which is endowed with the BI 
Lie bracket [X, F]n = XNF — yNX. The isomorphism property is a direct 
verification identical to Proposition 2.5 in |BlIsMaRa2005| . In particular, let 
{S, A, B), {S\ A', B') e (Sym(2X) (S) y) ® M and compute directly that 

^{[{S,A,B),{S\A\B')]) = 

= * (S" NS" - 5" NS', S NA' - S" N^, 1/4 ( M NA' - *A' NA)) 
' SNS' 





S'NS 




SNA' - 


-S'NA' 












S' NA) 


*ANA' - 


- ^A'NA 












N 0' 




'S' A'' 




'S' 


A'' 




'N 6 




's 


A' 









*A' AB' 




*A' 


AB' 











AB 



S A 
^A AB 

= [fiS,A,B),^iS',A',B')l^ 

as required. ■ 

Thus, we conclude that the BI system and the geodesic moment equations 
are equivalent. 



5.1 Klimontovich solutions of the Bloch-Iserles system 

The Klimontovich map for the geodesic Vlasov equation provides simple so- 
lutions of the Bloch-Iserles system in any dimension. For example, in even 
dimensions, the dynamics of the BI solution 

X{t) = Y.'^aC{t) (71) 

a 

is given by the system 

Ca=^'aJ2^bNClCa (72) 
b 

where N = 2J and ClCa is a covector such that (CtCa)j = iCDij iCaV ■ 
explicit index notation, one has 

(C)^-2^„^«;,rMCb),fe (O". 

b 

This system is a Hamiltonian system with the homogeneous quartic Hamiltonian 

1 

H^ = -Y.WaWtTv{\(:l)Cl) (73) 

a,b=l 



25 



Notice that, by writing the equation for (^^ as 

C, =WaJ2wtN(:lCa+ < lICa II' N (74) 

we can speciaUze the above to the simple case when = 1 for a fixed a and 
iffc = V6 a, so that 

x{t)^e{t) with C=IICII'nc (75) 

This case leads however to trivially linear dynamics, since the norm ||^|| is 
evidently conserved. This docs not happen for different norms in the quadratic 
moment Hamiltonian, such as H = 1/2 {y22 —I X'^ ® X^) ■ 

The above arguments also provide solutions to the Bloch-Iscrles system in 
2K + 1 dimensions. Indeed, the particle solution of EPCan (|55p becomes a 
solution of the BI system in the following form 

^(i) = E^' 

a 

Here Ca undergoes Hamiltonian dynamics with 

= i ^ u;a Wfc (C ■ Cfc + Tr ( \C) Cl) ) , (77) 

whose cohective EPCan equations are 

L=WaY.W,nClCa + \waY.W,nC,. (78) 

6 b 

As we shall see, these solutions are also momentum maps in both the even and 
odd-dimensional cases, since the Klimontovich solution (|55p is a momentum 
map. In particular, upon fixing A/" = 1, these are solution momentum maps 
V Sym(n) {C X), where is a symplectic space and X S Sym(n) is the 
Bloch-Iserles dynamical variable. This construction arises from a special case 
of the momentum map in ([55]) . with dim(S')=0 (Klimontovich case), V = V 
(symplectic vector space) and where the Lie algebra Xcan{V) is restricted to 
the Lie subalgebra of linear Hamiltonian vector fields. In the more general case 
when dim(S') > 1 one has also the conserved quantity in That is, the 

operation of taking moments preserves the dual pair structure of the Vlasov 
equation. 

Theorem 5.3 Upon fixing TV = 1, the solution ^76\ l of the odd- dimensional 
Bloch-Iserles system is a momentum map 

32K+1 ■■ {V,wl) ^Sym{2K + l) 

where is the symplectic form on the vector space V . Moreover, the solution 
|y_?[ ) in the even-dimensional case is also a momentum map 

i2K : {V,wJ) ^ Sym{2K) . 



Ca(i) 

4 



€ Sym(2is:+ 1) . 



(76) 
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Proof. In what follows we shall use the isomorphisms 



Sym(2A" + 1) ~ (Sym(2X) ® y) R 



and 



Sym(2i^) ~sp*(2ii',M), 



where sp{2K,M.) denotes the Lie algebra of Hamiltonian matrices. We prove 
the first statement, which comprises the second as a particular case. Let V be 
endowed with the Poisson structure 



then, the definition of momentum map can be verified by inserting G ~ (J,/3), 
with J = w;(C^C,l) e (Sym(2A')(S)F) ©M and (3 = (/Ja, /3i, /?o) its dual. Then, 



which identifies the infinitesimal action of linear (inhomogeneous) Hamiltonian 
vector fields on the phase space functions F e ^(V). 

Restricting to even 2K dimensions requires setting Pq = = f3i, thereby 
producing the action of (homogeneous) Hamiltonian vector fields, i.e. Hamilto- 
nian matrices in sp{2K,'R.). ■ 

The integrability properties of these solutions will be discussed elsewhere. 

6 Conclusions and outlook 

After reviewing the geometric basis of Vlasov moment dynamics, this paper 
showed how moment closures of the geodesic Vlasov equation (j26|) produce in- 
teresting known integrable systems, which include CH, CH2 and Bloch-Iserles 
(BI) equations in both odd and even dimensions. While the CH and CH2 
cases were already known to arise in ID jGiHoTr2005[ IGiHoTr2007] . the higher 
dimensional moment bracket showed that these moment closures also recover 
the EPDiff equation and extend CH2 to higher dimensions. The paper also 
recovered the BI system from (|26|) by a finite-dimensional moment closure, cor- 
responding to inhomogeneous quadratic phase-space functions. Thus, a kinetic 
theory approach led to special solutions of BI. 

The moment closures preserve the two equivariant momentum maps in the 
Vlasov dual pair. Preservation of this structure guaranteed that the resulting 
closed moment systems discussed here were still Poisson. This preservation also 
enabled reduction to finite-dimensional systems by using the Klimontovich par- 
ticle solutions from plasma theory. For example, the peakon solutions of the 
CH equation arose from a Klimontovich approach in |GiHoTr2005j . Singular 
solutions also arose upon allowing extra smoothing in the CH2 moment Hamil- 
tonian |GiHoTr2007[ IHoOnTr2009j . In addition, this paper showed in Section 
15.11 that the same approach also produced solutions of the finite-dimensional 




we have 
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BI system. However, Klimontovich solutions are not admitted by arbitrary ap- 
proximations. They are prevented, for example, when moment hierarchies are 
simply truncated at a certain weight. That is, moment closures preserve the 
Vlasov dual pair, while moment truncations do not, even though they may be 
shown to still be Lie-Poisson. Open questions concern both the construction 
of a Lax pair for the Klimontovich dynamics of the BI system and potential 
intcgrability properties of the truncated equations (e.g. the {Ai, A2) truncation 
for kinetic moments). 

The geometric setting showed how the left momentum map in the Vlasov 
dual pair recovers both the Klimontovich solution and the Lagrange-to-Euler 
map. That this geometry applies also to the Liouville equation illuminates the 
geometric footing of kinetic theory. Indeed, this paper explained how all the 
standard moment approximations in kinetic theory are momentum maps pre- 
serving the same dual pair. This construction would certainly be destroyed by 
introducing the collision integral, whose celebrated Boltzmann version implies 
irreversibility and produces a preferred direction of time via the 7J-theorem. 
This irreversibility prevents the Klimontovich solutions, which are solutions of 
a time-reversal invariant system. 

The kinetic-theory interpretation of geodesic Vlasov moment dynamics also 
provided insight into the physical description of the integrable cases. For ex- 
ample, the CH2 case was interpreted in this light as a charged fluid in the 
context of the one-component CH equation, thereby extending the CH model 
to include space-charge effects. On the other hand, CII2 has also been re- 
lated to shallow water dynamics by applying a series of approximations to the 
Green-Naghdi equations [Colv2008] . Remarkably, a modified version of CH2 
dynamics has also been applied in image matching [HoTrYo2007j . (In image 
matching, the Hamiltonian is the norm in which one applies optimal control.) 
Emergent pcakon solutions were found to result from applying smoothing 
to the Hamiltonian in |HoOnTr2009] . 

The paper also identified several other potentially interesting open problems. 
One of these is the problem of making physical applications of the Kirillov ad*- 
action for the Lie-Poisson bracket on the symmetric Schouten algebra for 
arbitrary values of {n,k). Another is to determine whether the Klimontovich 
average in plasma kinetic theory is a momentum map. One may also ask how 
the family of symplectically conserved quantities corresponding to statistical- 
moment versions of the Poincare invariants found in |HoLyScl990| may fit into 
the theory of kinetic moments. 

The kinetic approach used here may also provide physical interpretations of 
use in applying the BI system. For example, particle beams in linear accelera- 
tor lattices are described in terms of symplectic transfer matrices. (The same 
holds for linear ray optics.) In formulas, one has the relation z{t) = M{t) z(0), 
where A4{t) is a one parameter subgroup of Sp(6,K) determining the beam 
evolution z{t) C K^. In this sense, geodesies in Sp(6,]R) would correspond to 
optimal transfer maps for particle or optical beams. It is interesting that similar 
approaches have recently emerged in quantum computation, where the Hamilto- 
nian of the system is constrained to an optimal trajectory (i.e., a geodesic) by a 
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cost function from optimal control theory [BrElHoOSl INiDoGuDo08| . These ad- 
ditional open problems bode well for the potential success in future applications 
of using the geometric approach to Vlasov dynamics discussed here. The appli- 
cation of ideas from optimal control to the Vlasov moments may be especially 
fruitful. 
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